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Abstract 

The scalar strange-quark matrix element of the nucleon is computed with lattice QCD. A mixed- 
action scheme is used with domain-wall valence fermions computed on the staggered MILC sea- 
quark configurations. The matrix element is determined by making use of the Feynman-Hellmann 
Theorem which relates this strange matrix element to the change in the nucleon mass with respect 
to the strange quark mass. The final result of this calculation is m s (N\ss\ N) = 49 ± 10 ± 15 MeV 
which can be translated into f s = m s (N\ss\N) /m^ = 0.053 ± 0.011 ± 0.016. 

Given the lack of a quantitative comparison of this phenomenologically important quantity 
determined from various lattice QCD calculations, we take the opportunity to present such an 
average. The resulting conservative determination is (all uncertainties added in quadrature) 
f 8 = 0.043 ± 0.011 (99% confidence interval). 
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I. INTRODUCTION 



Determining the strange content of the nucleon has been a long standing interest of nuclear 
and particle physicists. The scalar strange content of the nucleon can be related to kaon- 
nucleon scattering and phenomenologically to the possible condensation of kaons in dense 
nuclear environements [U [2] . The strange content of the nucleon may also play an important 
role in the scattering of dark-matter particles off nuclei. The general low-energy coupling of 
dark matter to nuclei has recently been worked out systematically using low-energy effective 
field theory [31 H] . The spin-independent coupling is the simplest and hence received the most 
attention historically. The spin-independent elastic scattering of dark-matter off a nucleon 
is proportional to the square of the scalar matrix elements (N\m q qq\N) for quarks of flavor 
q. [SHU]. There are no direct experimental means of measuring these matrix elements. The 
heavy quark q = {c, b, t} matrix elements can be computed from perturbative QCD and are 
reasonably well known [10J. The light quark q = {u, d} matrix elements can be reasonably 
determined from nN scattering [H] [12]. The scalar strange quark matrix element presents 
the most theoretical challenge to determine reliably and has contributed one of the largest 
uncertainties in dark-matter detection experiments [6]. For these reasons, there has been a 
resurgent interest in determining m s (N\ss\N) using lattice QCD, beginning with the work 
in Refs. [U [13] . It is more common in the context of dark-matter searches to normalize this 
quantity by the nucleon mass, 

m 8 (N\ss\N) 
m N 

There are two typical approaches used to determine this quantity from lattice QCD. The 
scalar strange quark matrix element can be directly computed or one can take advantage of 
the Feynman-Hellmann Theorem; 

m s (N\ss\N) = . (2) 

Most groups use the direct method [TH - 120] . one group uses a hybrid approach which involves 
elements of both methods [2"T| |2"2"] , and some groups use the Feynman-Hellmann method [TBI 
Unj |2"3T[2"T] . For a recent review of the scalar strange content of the nucleon, see Ref . [2E] • 

The present work utilizes the Feynman-Hellmann Theorem which has some distinct ad- 
vantages over the direct method: it is numerically less expensive and the ground state contri- 
butions to the two-point correlation functions can be significantly more reliably determined 
than plateaus in direct matrix element calculations. 

We begin by presenting details of our lattice calculation in Sec. [IT] and then present 



the determination of m s (N\ss\N) in Sec. Ill We have found a quantitative comparison 



of various lattice QCD calculations of this quantity lacking in the literature. Given its 
important phenomenological role, we were compelled to compile such a comparison which 
we provide in Sec. |IV| along with the results of the present work. While a precise and reliable 
determination of f s is still lacking, there is a welcoming consistency in the determination of 
this quantity from a wide variety of lattice calculations. 

II. DETAILS OF THE LATTICE CALCULATION AND NUMERICAL RESULTS 

The present work utilizes mixed-action lattice QCD calculations with domain-wall 
fermion [29H33] propagators computed on the rif = 2 + 1 asqtad-improved [MI 35J rooted, 
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staggered sea-quark configurations generated by the MILC Collaboration [361 E7J, (with 
hypercubic-smeared [38H1T] gauge links to improve the chiral symmetry properties of the 
domain- wall propagators), a strategy initiated by the LHP Collaboration [I2H37]. A prin- 
ciple motivation for this choice is the good chiral symmetry properties of the domain-wall 
action, while utilizing the less numerically expensive lattice configurations of the staggered 
action. It has been shown that the chiral symmetry properties of the valence domain- 
wall fermions highly suppresses sources of chiral symmetry breaking from the sea-quark 
action [^8H5Tj . This has led to a number of important results including a determination 
of the kaon bag parameter Bk [52] , the charmed and static baryon spectrum [531 El] , hy- 
peron axial charges [55], as well as a number of results from the NPLQCD Collaboration 
including two-hadron scattering lengths [561462] . multi-meson interactions, condensates and 
the three-pion interaction [5"3TI5"5] . as well as a number of others [66-69J. There have been a 
few other choices for mixed actions all utilizing overlap [TO] ITT] valence-fermions on a vari- 
ety of sea-quark configurations: Wilson sea-fermions [72] , twisted-mass sea-fermions [73] [71] , 
domain- wall sea fermions [75j[76], and HISQ sea fermions [77]. Mixed-action calculations are 
inherently unitarity violating with partially quenched effects only vanishing in the continuum 
limit. It is therefore imperative to compare numerical results with the scaling violations pre- 
dicted from the mixed-action effective field theory [18H5"T| [751487] . This has been undertaken 
in a exploratory extent with baryons [HI [88] but the only systematic studies have been with 
the do correlator [BUI [HTJ, which is highly contaminated by the unitarity violating effects, 
and a recent determination of low-energy-constants in the two-flavor chiral Lagrangian for 
pions [BH]- Despite the limited study of discretization effects, there are reasons to believe 
they are small for many quantities [l6j H8H50[ 190] . 

A. Parameters of the Lattice QCD Calculation 

The present calculation utilizes the Feynman-Hellmann Theorem to determine the scalar 
strange quark matrix element in the nucleon, Eq. ([2]), limiting the work to a small set of 
available ensembles. Details of the various ensembles and parameters are collected in Table [1} 
There are two sets of ensembles at the h « 0.125 fm lattice spacing with fixed light quark mass 
and strange quark masses that straddle the physical strange quark mass. These are denoted 
by the sets m sea = {m010m030, m010m050} and m sea = {m030m030, m030m050} respec- 
tively. On the b « 0.09 fm ensembles, there are two sets with fixed light quark mass and 
strange quark masses straddling the physical strange quark mass. In this work, preliminary 
results are presented only for one of these sets with m sea = {m0031m0186, m0031m031}. 
The values of the domain-wall quark masses, the fifth dimensional extent L 5 and the domain- 
wall mass M5 were taken from the NPLQCD production runs [89] . 

B. Results of the lattice calculation 

The light and strange quark propagators were constructed with a Gaussian-smeared 
source [9U E2] and both smeared (SS) and point (SP) sinks. Correlation functions were 
then constructed with the quantum numbers of the pion and proton. The pion masses were 
determined with a fully correlated simultaneous fit to the SS and SP correlation functions, 
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TABLE I: Lattice parameters used in the present work. For some of the calculations, the time 
direction was chopped at t = 32 with Dirichlet boundary conditions (denoted by volumes with 
x32). 











bz 


■i 0.125 fm ensembles 








p 




V 


M 5 


L 5 


i dw f 

bm^ 


6m[ es 


i dw f 

bm s 


bm r s es 


N src X N cfg 


6.75 


m010m030 


20 3 x 64 


1.7 


16 


0.0138 


0.001564(03) 


0.081 


0.000892(2) 


53 x 328 


6.76 


m010m050 


20 3 x 64 


1.7 


16 


0.0138 


0.001566(11) 


0.081 


0.000913(2) 


4 x 656 


6.76 


m010m050 


20 3 x 32 


1.7 


16 


0.0138 


0.001552(27) 


0.081 


0.000913(2) 


24 x 769 


6.79 


m030m030 


20 3 x 64 


1.7 


16 


0.0478 


0.001052(04) 


0.081 


0.000809(4) 


30 x 367 


6.81 


m030m050 


20 3 x 32 


1.7 


16 


0.0478 


0.001013(06) 


0.081 


0.000862(7) 


24 x 564 










b'. 


« 0.09 fm ensembles 








P 


f^sea 


V 


M 5 


L 5 


, dw f 

bm l 


bm\ es 


bm s 


bm r s es 


N src x N cfg 


7.08 


m0031m0186 


40 3 x 96 


1.5 


12 


0.0035 


0.000431(3) 


0.0423 


0.000236(2) 


1 x 356 


7.08 


m0031m031 


40 3 x 96 


1.5 


12 


0.0035 


0.000428(3) 


0.0423 


0.000233(2) 


1 x 422 



with a single cosh used for both correlators; 

C (S x)(t) ~ A (sx) e- m * T ' 2 cosh(m^ - T/2)) , (3) 

where X = S,P. In all cases, the configurations are binned until the statistical uncertainty 
stoped changing appreciably. To determine the fitting systematic, the length of the time- 
extent used in the fit and the starting time were varied over a wide range, with a minimum 
plateau length of ~ 0.5 fm. For each fit, the Q- value was used as a weight, where 

roc 

Q = / d X 2 V( X 2 ,d), (4) 
J-x 2 

with the probability distribution function for \ 2 with d degrees of freedom 

The results of these fits are plotted over a representative window in time along with cosh- 
style effective masses 



cosh/ 
eff 



, 1 , i fC(t + r) + C(t-T)\ 

in Figures |l| and H The (black) squares are from the SP correlation functions while the 
(colored) open circles are from the SS correlation functions. The right-side bar in each plot 
displays the mass probability distribution function determined from 



To determine the central value we take 



TTliQi , . 
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FIG. 1: Pion mass effective mass plots on the b ~ 0.125 fm ensembles. 



In many cases, the systematic is approximately Gaussian, and so the 16% and 84% quantiles 
are used to determine the systematic uncertainties. In all plots, the inner (colored) band 
represents the statistical uncertainty in the fit while the outer band represents the statistical 
and fitting systematic added in quadrature. In the case of the m sea = m030m050 ensemble, 
only results with Dirichlet boundary conditions in time are available. For this case, the 
correlation functions are fit simultaneously with a single exponential 

C(sx)(t) ~ A( S x) exp(-m 7r t) , (9) 
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FIG. 2: Pion mass effective mass plots on the b « 0.09 fm ensembles. 



TABLE II: Computed pion and nucleon 
of T\jb used to convert to physical units 



masses on the various ensembles. Additionally, the value 
is provided, obtained from Refs. [37^ [93], 



P 



m. 



V 



brriN 



r f(bmf y ,bmf\(3) 



6.75 


m010m030 


20 3 x 64 x 


16 


0.22178(33) (Ij) 


0.7177(18)(i§) 


2.711(4) 


6.76 


m010m050 


20 3 x 64 x 


16 


0.22285(28) (f) 




2.739(3) 


6.76 


m010m050 


20 3 x 32 x 


16 




0.7311(19)(i) 


2.739(3) 


6.79 


m030m030 


20 3 x 64 x 


16 


0.37323(27)(20) 


0.8653(17)(i) 


2.821(7) 


6.81 


m030m050 


20 3 x 32 x 


16 


0.37493(26) Qf) 


0.8740(18)(H) 


2.877(4) 


7.06 


m0031m0186 


40 3 x 96 x 


12 


0.10192(38)(1) 


0.4621(64)(|§) 


3.687(4) 


7.08 


m0031m031 


40 3 x 96 x 


12 


0.10165(35)(f|) 


0.4603(48) (Ijl) 


3.755(4) 



and compared with the standard effective mass 

1 



,ln 
eff 



{t > T) = T l \C(t + T) 



C(t) 



(10) 



Clearly, in this case, the ability to explore the fitting systematic is more limited. For reasons 
discussed in Ref . [89] , the dip in the effective mass is not believed to represent a lower ground 
state energy, but rather contaminations from the Dirichlet boundary condition. The results 
are collected in Table HT1 

The proton masses are trickier to determine as the signal-to-noise ratio decays exponen- 
tially in time [94^ 



lim 

t^oo N(t) 



(11) 



The mass determined in a given fit from tj to tf is then susceptible to larger fitting sys- 
tematics. It is not uncommon for the effective mass plateau to shift by order one standard 
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deviation and form a new plateau either higher or lower at times when the statistical fluctu- 
ations grow appreciably. It is therefore important to develop a systematic analysis algorithm 
that both takes advantage of the precise statistical fluctuations at early times while allowing 
for the possibility that the late time fluctuations represent the true groundstate. It was 
demonstrated that correlation functions determined with (9(10 5 ) reasonably statistically in- 
dependent sources on 0(1O 4 ) Monte-Carlo trajectories, a variety of analysis methods could 
be used all producing consistent results [95]. With fewer measurements, not all methods 
work as well. 

One technique which works better than others is the Matrix-Prony method, as described 
in Refs. [951 96] . With this method, the ansatz is made that the correlation functions can 
be described with a transfer matrix 

y(t + r) = f(r)y(t), (12) 

where in our case 

One solves for T(r) = M~ l (r)V by constructing the matrices 

/ t +At \ _1 /t +At \ _1 

V t=to / V t=t J 

One then solves the eigenvalue equation for the principle correlators 

T(r)q n = {\ n ) T q n , with A„ = e~ E - . (15) 

The original ansatz is satisfied if over the range of time to to to + At, the resulting principle 
correlation functions are well described by a single exponential. 

In this work, to determine the fitting systematic, the choices of to and At are varied over 
a wide range, with At > 0.5 fm. For each choice, the ground state principle correlation 
function is fit with a single exponential, Eq. (j9|, over a range of time t{ — tf, chosen indepen- 
dently of to and At. The initial and final times in the fit are also varied over a wide range 
under the constraint tf — ti > 0.5 fm. For each fit, the Q- value is recorded along with the 
statistical uncertainty of the fit. The various fits are then averaged with the weight similar 
to that of the pions, but also suppressed by the statistical uncertainty of the fit; 

m = — ^ with Wi = — . (lb) 

In this way, the plateaus at later times, with larger uncertainties, and hence larger Q- values, 
do not dominate the determination of the fitting systematic. The resulting fits are displayed 
along with effective mass plots of representative Matrix-Prony determinations of the ground 
state principle correlation function in Figures [3] and |4j In these figures, the colored effective 
mass points represent the time window over which the Matrix-Prony method is applied while 
the gray effective mass points lie outside this region. As is evident, the resulting systematic 
mass-probability distribution tends not to be Gaussian. For simplicity, we still take the 16% 
and 84% quantiles to define the systematic uncertainty. The inner colored bands represent 
the statistical uncertainty and the outer gray bands represent the statistical and systematic 
uncertainties added in quadrature. 
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FIG. 3: Proton mass and representative effective mass plots on the b ~ 0.125 fm ensembles. 



C. Scale setting 



To convert these lattice quantities to physical units we use the scale setting procedure 
described in Ref. [89]: the dimensionless lattice results are converted into r\ units, with 
y (bmi, bm s , (5) determined by the MILC Collaboration on each ensemble, and extrapolated 
to the physical values of the light quark masses [33 E3] , which are listed in Table [XTJ The 
MILC Collaboration has also determined the physical value of r\ [371 E3] , 



r x = 0.31174(20) fm. 



(17) 



which is used to then convert all values into physical units, Table III 
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FIG. 4: Proton mass and representative effective mass plots on the b ~ 0.09 fm ensembles. 



TABLE III: Computed masses and decay constants converted to MeV with r\ = 0.31174(20) fm. 



p 




V 


[MeV] 


m N [MeV] 


6.75 


m010m030 


20 3 x 64 x 16 


380.5(.6)(;|) 


1231(3)( 3 ) 


6.76 


m010m050 


20 3 x 64 x 16 


386.3(.5)(:|) 




6.76 


m010m050 


20 3 x 32 x 16 




1267(3)(«) 


6.79 


m030m030 


20 3 x 64 x 16 


666.4(.5)(.4) 


1545(3)(1) 


6.81 


m030m050 


20 3 x 32 x 16 


682.7(.5)(;i) 


1591(3)( 6 7 ) 


7.06 


m0031m0186 


40 3 x 96 x 12 


237.8(0.9)(1.3) 


1078(15)(||) 


7.08 


m0031m031 


40 3 x 96 x 12 


241.6(0.8) (?;°) 


1094(11)G|) 



There is an important additional advantage to our method of scale setting. To invoke 
the Feynman-Hellmann Theorem, the change in the nucleon mass with respect to a change 
in the strange quark mass must be undertaken with all other parameters held fixed jZQ [22] . 
By converting to physical units with values of ri/b extrapolated to the physical light and 
strange quark masses, this is achieved with only a small mistuning of the light quark mass. 
This mistuning can be determined by multiplying the values of r\jb in Table |LT| with the bare 
input light quark masses, which is found to be < 1%. The difference in the pion masses in 



physical units, Table III (for the same light quark masses) is due in part from the different 
strange quark masses and also in part from the light-quark mass mistuning. Using leading- 
order chiral perturbation theory (%PT), m\ = 2fi ^i«, the difference in the pion masses 
is larger than expected from just the mistuning of bmi, but of comparable magnitude. A 



simple average of the pion masses in Table III is used. 
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TABLE IV: Extracted values of mjv(mf y ) and mf v (N\ss\ N). The first uncertainty is statistical, 
the second fitting systematics and the third is from the uncertainty on the determination of m p s hy . 



m n [MeV] 


m N {m phy ) [MeV] 


ml hy \N\ss\N) [MeV] 


383.4(.6)(;§) 


1243(2)(3)(1) 


63(8)(11)(1) 


674.6(.5)(.4) 


1560(2)(4)(2) 


82(8)(14)(2) 


240(1)(2) 


1107(21)(33)(1) 


72(54)(89)(1) 



III. THE STRANGE SCALAR MATRIX ELEMENT IN THE NUCLEON 



As discussed in the introduction, there are two methods for determining the scalar strange 
quark matrix element in the nucleon: a direct calculation of the matrix element employed 
by some groups [T4T - [20] and an indirect determination through the Feynman-Hellmann The- 
orem [131 1231427] . Eq. (J2|, and a hybrid approach [211 E2]. This work utilizes the latter 
method. For each light quark mass ensemble, we have a determination of the nucleon mass 
at values of the strange quark mass which straddle the physical strange quark mass. These 



results, Table III, can be used to interpolate to the physical value of the strange quark mass, 
Taylor expanding about bm p s hy ', and determine the two quantities 



m N (m phy ), 



dm N (m s ) 



dm* 



(18) 



The MILC Collaboration has determined values of the strange quark mass to be bm phy = 
0.0350(7) and bm phy = 0.0261(5) on the b m 0.125 fm and b m 0.09 fm ensembles respec- 
tively [371 197]. Performing the Taylor expansion about m phy , the values of the nucleon mass 



and m s (N\ss\N) are determined on each ensemble. These results are collected in Table IV 
and the resulting interpolations are displayed in Figure [5j In these figures, the vertical 
dashed lines represent the 68% confidence interval for the determination of bm phy on the 
b ~ 0.125 fm and b ~ 0.09 fm ensembles. The uncertainty on bm phy is included in the anal- 
ysis and represented by the third uncertainty in Table IV On the b « 0.125 fm ensembles, a 
precise determination of the scalar matrix element is obtained. However, on the b ~ 0.09 fm 
ensembles, the results are too imprecise to determine a non-zero value. 



A. Chiral Extrapolation 

The results for m s < A^|ss| N > must be extrapolated to the physical value of the pion mass. 
In Ref. [HH], the two flavor extrapolation formula for this matrix element was determined at 
next-to-leading order (NLO) in the chiral expansion, 

(N\ss\N) =< N\ss\N >° -^p (< N\ss\N >° - < A\ss\A >°) j£ + E s ^~ 2 , (19) 

where (H\ss\H)° represent the leading order (LO) contribution to the scalar strange matrix 
element in the hadron H, g n NA is the axial pion-nucleon-delta coupling appearing in the 
SU(2) baryon chiral Lagrangian, is a chiral loop function non-analytic in the pion mass 
and the delta-nucleon mass splitting (A = wa — mjv) and E s is a low energy constant 
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FIG. 5: Nucleon mass versus the strange quark mass on the b ~ 0.125 fm and b ~ 0.09 fm 
ensembles. The vertical dashed lines represent the 68% confidence interval for the determination 
of bmP s hy on the b ~ 0.125 fm and b ~ 0.09 fm ensembles. 



appearing at NLO. In the large- N c expansion, the LO matrix elements for the nucleon and 
the delta are both O(N^) but there is no cancellation at this order [99], so one does not 
expect a strong cancellation between these NLO contributions. 1 In principle, one should 
use the partially quenched formula, also provided in Ref . [98j , and convert it to the relevant 
mixed-action formula [50] to perform the extrapolation. However, clearly the most significant 
shortcoming of the present work is the limited number of light quark mass points. With 
non-zero results at only a single lattice spacing, the mixed-action extrapolation can not 
be performed regardless. The best that can be done with the present results is a simple, 
effectively zero-degree of freedom extrapolation using the formula 

m s (N\ss\N) = c + c 2 ml . (20) 

While this will not result in a precise and accurate determination of the scalar strange matrix 
element, it will provide a good guide to the approximate value at the physical point. While 
not a rigorous expectation, it has been found that matrix elements of the nucleon tend to 
have very mild pion mass dependence, see for example the recent review [101] . Performing 



See also Ref. [100 for further discussion on the baryon massess in the large iV c counting. 
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FIG. 6: Light quark extrapolation of ms(N\ss\N) versus m 2 . The location of the vertical dashed 
line is given by (m£ y ) 2 . 



this simplistic pion mass extrapolation, using the isospin averaged y = 138.0 MeV, 

m p s hy (N\ss\N) h = 56 ± 12 ± 17 MeV. (21) 

The extrapolation is displayed in Figure |6j 

Given the limited ability to perform the chiral extrapolation, we also explore the light 
quark mass dependence of f s = m s {N\ss\N) /m^ to improve the estimate of systematic 
uncertainties. It has been observed that the nucleon mass displays a remarkably linear 
dependence on the pion mass [HIl EO] • For this reason, two extrapolation functions are used 
to estimate extrapolation systematics: 

/- = /i 0) + fPml , (22a) 

/. = /i 0) + /i 1 V , (22b) 

yielding the results 

f a = 0.050 ±0.009 ±0.014, (23a) 

f 8 = 0.050 ±0.012 ±0.018, (23b) 

respectively. These extrapolations are displayed in Figure [7j The quantity f s is observed to 
have negligible light quark mass dependence. 

These results can be compared with the extrapolation of m s (N\ss\N) by converting 
with the isospin averaged nucleon mass m p ^ y = 938.9 MeV. In Table 



V 



these three different 

extrapolation results are collected. Additionally, a correlated weighted average is performed. 
To perform the correlated average, Gaussian distributions of the results in Table IV are 



created independently for each light quark mass point, with Ncauss = 10 4 in all cases. For 
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FIG. 7: Extrapolation of f s . The location of the vertical dashed line in each plot is determined 
from m^ hy . 



TABLE V: Extrapolated values of m s {N\ss\N) and f s . These results are averaged in a weighted 
and correlated fashion described in the text. 



Quantity Extrapolated Extrapolation Function m% y (N\ss\ N) [MeV] 



m p s hy (N\ss\N) 

fs 
fs 



Eq. (20) 



Eq. (22a) 



Eq. (22b) 



56 ± 12 ± 17 
47 ± 9 ±13 
47 ± 12 ± 17 



0.059 ±0.012 ±0.019 
0.050 ±0.009 ±0.014 
0.050 ±0.012 ±0.018 



Correlated Average 



49 ± 10 ± 15 



0.053 ±0.011 ±0.016 



each sample, all three extrapolations are performed, preserving the correlations between 
the fits, with inverse weights given by the statistical and systematic uncertainties on the 
individual mass points. For each sample, these three results are then averaged with weights 
given by the inverse uncertainties from the individual analyses (quoted in Table \V\ . This 
yields the final result 

m s (N\ss\N) =49 ±10 ±15 MeV, (24a) 
f 8 = 0.053 ±0.011 ±0.016. (24b) 



IV. RESULTS AND DISCUSSION 



For the present work, the Feynman-Hellmann Theorem was invoked to determine the strange 
content of the nucleon through a change as the strange quark mass is varied 

/An- I ,r\ dm N 
m s (N\ss\N) = m s-g^- ■ 

By taking care to set the scale using values of r\jb which were extrapolated to the physical 
values of the light and strange quark masses, the nucleon mass variation was determined 
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with all other parameters held constant (with precision better than 1%), as is required 
for a proper determination of this quantity [211 122] • There are several groups who have 
used the Feynman-Hellmann Theorem [T3"l IT9"1 I23H2Z] as well as more determinations with 
a direct calculation of the matrix element [T4T[2"0"] and results from a hybrid approach [2"Tl 
[22] • Before making a detailed comparison with other works, we first highlight advantages 
and disadvantages of the present work. The distinct advantage of using the Feynman- 
Hellmann Theorem over direct methods is that the ground state plateau of the nucleon can 
be significantly more reliably determined than the plateau for the matrix element calculation, 
see the plots of ratio determinations in any of Refs. [T4T420"] (the direct calculation requires a 
vacuum subtraction, adding substantial statistical noise). The disadvantage of most groups 
employing the Feynman-Hellmann Theorem is the reliance upon SU (3) baryon %PT [TBI I25T 
|2"T] , which is known to not have a converging expansion for the nucleon mass |4"o] E21 EH1 11021 
103J. Therefore, it is not clear the full extrapolation systematic has been properly addressed 
in those works. 2 This concern is substantiated by the discrepancy between independent 
SU(3) baryon x;PT analyses and their determination of f s [l_3j I25T - T2T] . 3 The current work 
does not suffer from this issue. 

The most severe limitation of the present work is the small number of light quark mass 
points (two) for which there is a non-zero determination of m s (N\ss\N) . Given the sig- 
nificant numerical cost of the domain-wall propagators on the b rj 0.09 fm ensemble with 
m v ~ 240 MeV, it is not clear how soon a more precise determination will be obtained at 
this point. Given the very mild light quark mass dependence observed in this work, and in 
nucleon matrix elements in general, we believe the present determination offers a reliable 
estimate of the scalar strange content of the nucleon, but neither a precise nor demonstrably 
accurate value. Our final result is 

m s (N\ss\N) =49 ±10 ±15 MeV, 

f s = 0.053 ±0.011 ±0.016. 

A. Lattice QCD Comparison and Average 

Given the phenomenological importance of the scalar strange content of the nucleon, see 
for example Refs. [5H9], it is prudent to review the limitations of the present determination 
and to compare and contrast these results to other lattice QCD determinations. There are 
two results which use the same MILC ensembles with staggered valence quarks [2~Tl [22] and 
one determination with the same mixed-action scheme but a direct determination [20]- It is 



2 The work in Ref. [53] also uses SU(3) baryon %PT, but uses a variety of other extrapolation methods, 
resulting in a conservative estimate of their uncertainties. 

3 Despite these criticisms, we point out in Ref. [26;, a striking agreement is found between baryon mass 
results extrapolated from one set of lattice calculations 1102] . with SU(3) baryon %PT, and then used 
to predict results from a completely independent calculation [lOlj. Moreover, independent verification 
of the consistency of various lattice calculations of the ground state baryon spectrum and SU(3) baryon 
%PT has been found independently [251 11051 1106j . 
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_A 

0.063(11) 
0.032(25) 
0.012O 
0.014(06) 
0.048(15) 



[21]n/ = 2 + l 
[14] n, = 2 
[16] n, = 2 
[17]n/ = 2 + l + l 
[18] n/ = 2 + 1 



CD 
Q 



0.009(22) 
0.046(11) _ 
"0.058(09) 



[19]n/ = 2 + l 
[20] n f = 2 + l 



[22]ra/_=2 + l 

'[19]n/ = 2 + l 
[13] ra/ = 2 + 1,5(7(3) 
[23] ra/ = 2 + 1 
[24]n / = 2 + l 
[25] n,= 2 + 1,5(7(3) 
[26] r( / = 2 + 1, 5(7(3) 
[27] -;(/ = 2 + 1,5(7(3) 
present work 



J 

"55 
X 
c 

£ 



0.023(40) 
0.033(17) 
0.036(tg) 
-0,075(73) 
0.023(22) 
0.022(lf) 
0.131(63) 



0.053(19) 



0.043(11) 



lattice average (see text) 



(1.00 



(l.0 r ) 



0.10 



FIG. 8: Comparison and average of lattice QCD calculations of f s as described in the text. The 
quoted uncertainties are taken as the statistical and systematic uncertainties added in quadrature 
from a given reference, nj = 2 + 1 indicates a dynamical strange quarks as well as up and down. 
SU(3) is used to indicate results which rely heavily on SU{3) baryon xPT. Some results are 
excluded for various reasons but displayed to demonstrate their consistency: [21] was updated in 
[22] , the rif = 2 results were not averaged with the nj = 2 + 1 [TJ1 [16] , the results in [17] were 
preliminary and not extrapolated to the physical pion mass, the results in [18] are preliminary and 
only exists in a conference proceedings. 

interesting to first compare our results with these; 



m s (N\ss\N)[Me~V] 



59 ± 6 ± 8 Ref. (21 
54 ± 5 ± 6 Ref. [22 



I3±S + G Ref. 120 ' ^ 



49 ± 10 ± 15 present work 



In the literature, there is currently no determination of f s which considers all the available 
results from lattice QCD, and so we take the opportunity to provide one here. 4 We use an 
approach similar to the FLAG working group of FLAVIANET which has provided lattice 
determinations of various quantities important to low-energy hadronic physics [107] . In 
particular, the FLAG working group has developed a scheme to judge the confidence to 
place in various determinations, based upon standards such as the lightest pion mass used, 
whether or not a continuum limit has been performed, and whether the infinite volume limit 
has been performed. For each criterion, a green star (~k) is awarded to results which meet 
the strictest constraints, and orange circle (•) is given to results with room for improvement 
and a red square (■) to those with room for significant improvement. This provides a useful 



There is a recent review on the topic in Ref. [55], but a lattice average is not provided. 
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guide to people outside the lattice community and motivation for those in the community 
to improve their results. 

Using the standards of Ref. |lU7j . most results for f s receive an orange circle. There is 
one group which receives the green star and the rest receive a red square. The results with 
a red square suffer either from too few light quark mass points to make a reliable chiral 
extrapolation or rely to heavily on SU(3) baryon x?T. There are two analyses which we 
promote from a red square to an orange circle because while they rely heavily on SU(3) 
baryon xPT, they have demonstrated a remarkable consistency of their analysis with four 
or more independent lattice calculations [251 126] . We exclude results which are either not 
published or not in an arXiv posting (as results in conference proceedings often undergo- 
ing larger-than-quoted systematic changes). We further exclude results which have yet to 
perform light-quark mass extrapolation and results calculated without dynamical strange 
quarks (rif = 2) are not included in the average. These results are displayed in Figure [8] 

For the scalar strange content of the nucleon, the current state of results is such that a 
simple weighted average of good (green star) results can not be performed in a meaningful 
way. As can be seen in Figure [HJ there is good consistency between most of the results. 
There are not a large number of orange circle results, so we chose to include all results in 
the average. Moreover, we believe despite their red-square assignment, these results offer 
valuable information which should not be ignored at this time. 

A simple weighted average, using the quoted uncertainties as the inverse weights, pro- 
duces an unbelievable small final uncertainty. This also ignores the fact that systematic 
uncertainties are typically non-Gaussian, and in the case of lattice QCD calculations, not 
cleanly separable from the statistical uncertainties. Moreover, it does not account for the 
quality of the results, judged using the rubric of the FLAG working group. In an attempt 
to include all these issues, the following ad hoc procedure is used to perform a weighted 
average of all the results (presented in Figure [8]); 

i) for each result, an independent random sample is generated with a sample size of 10 4 , 
drawn from a uniform distribution between the quoted uncertainties 

ii) for each random sample, a weighted average of all results is performed, with weight 

Wi = yi/ai , 

where <7j is the total (symmetric) uncertainty from a given result, and we arbitrarily 
chose Hi = 1,2,3 for the red square, orange circle and green star respectively: an extra 
reduction of 0.5 is assigned to results which rely heavily on SU(3) baryon %PT 

iii) the mean and 99% confidence intervals of the resulting distribution are quoted, see 
Figure [8] 

This yields a lattice average value of 

m s (N\ss\N) = 40 ± 10 MeV, 
f s = 0.043 ±0.011. 

As was first discussed in Refs. [El [13], there is now compelling evidence from lattice QCD 
that the value of the scalar strange content of the nucleon is substantially smaller than 
previously estimated and does not play as significant a role in dark-matter searches as 
previously thought [6j. For a recent review of the lattice QCD determinations of the scalar 
strange content of the nucleon, see Ref. [28] . 
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